Rank Restriction for the Variational Calculation of Two-electron Reduced Density 
Matrices of Many-electron Atoms and Molecules 
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Variational minimization of the ground-state energy as a function of the two-electron re- 
duced density matrix (2-RDM), constrained by necessary A'^-representability conditions, provides 
a polynomial-scaling approach to studying strongly correlated molecules without computing the 
many-electron wavefunction. Here we introduce a new route to enhancing necessary conditions for 
A''-representability through rank restriction of the 2-RDM. Rather than adding computationally 
more expensive A'^-representability conditions, we directly enhance the accuracy of two-particle (2- 
positivity) conditions through rank restriction, which removes degrees of freedom in the 2-RDM 
that are not sufficiently constrained. We select the rank of the particle-hole 2-RDM by deriving the 
ranks associated with model wavefunctions including both mean-field and antisymmetrized geminal 
power (AGP) wave functions. Because the 2-positivity conditions are exact for quantum systems 
with AGP ground states, the rank of the particle-hole 2-RDM from the AGP ansatz provides a mini- 
mum for its value in variational 2-RDM calculations of general quantum systems. To implement the 
rank-restricted conditions, we extend a first-order algorithm for large-scale semidefinite program- 
ming. The rank-restricted conditions significantly improve the accuracy of the energies; for example, 
the percentages of correlation energies recovered for HF, CO, and N2 improve from 115.2%, 121.7%, 
and 121.5% without rank restriction to 97.8%, 101.1%, and 100.0% with rank restriction. Similar 
results are found at both equilibrium and non-equilibrium geometries. While more accurate, the 
rank-restricted A''-representability conditions are less expensive computationally than the full-rank 
conditions. 

PACS numbers: 31.10.-(-z 



I. INTRODUCTION 

Significantly more information is encoded within the 
wavefunction than is necessary for the calculation of en- 
ergies and properties of many-electron quantum systems. 
In 1955 Mayer proposed in Physical Review calculating 
the ground-state energy variationally as a functional of 
the two-electron reduced density matrix (2-RDM) rather 
than the iV-electron wavefunction [ll-Q . Unlike the wave- 
function the 2-RDM scales polynomially with the num- 
ber N of electrons in the molecule. In further work, 
however, it became apparent that the 2-RDM must be 
constrained by non-trivial conditions to ensure that it is 
representable by an A'^-electron density matrix (or wave- 
function), and the search for these conditions became 
known as the N-representability problem [1, 0-^3] ■ For 
nearly 50 years the direct calculation of the 2-RDM with- 
out the wavefunction was stymied by the need for bet- 
ter 7V-representability conditions and better optimization 
methods. 

The variational computation of an A^-particle system's 
ground-state energy as a functional of the 2-RDM has 
recently been realized through advances in (J) develop- 
ing 7V-representability conditions d, H, H, |^ IT^ I22I 
and {ii) designing optimization algorithms |12h151 |25| - 
[28j . A systematic hierarchy of A^-representability condi- 
tions has been developed in the form of p-positivity con- 
ditions [sl. [qI. [ill [l7| ■ which constrain p -|- 1 distinct metric 
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matrices of the fermionic p-RDM to be positive semidef- 
inite (A matrix is positive semidefinite if and only if its 
eigenvalues are nonnegative) . The p + 1 metric matrices 
represent the probability distributions of p — g particles 
and q holes with q ranging from to p where a hole 
is the absence of a particle '9]. The p-positivity condi- 
tions ensure that each of these probability distributions 
is nonnegative. These conditions, even for small p, are 
capable of capturing both moderate and strong electron 
correlation; for example, the 2-positivity conditions are 
necessary and sufficient for computing the ground-state 
energies of pairing Hamiltonians JTlj, often employed in 
modeling long-range order and superconductivity. 

Minimizing the ground-state energy as a 2-RDM func- 
tional constrained by these conditions produces a special 
type of optimization known as semidefinite programming 
(SDP) [IMI,!!!-!!!. importantly, because SDP prob- 
lems are solvable in polynomial time, the variational 2- 
RDM method provides a polynomial-time relaxation of 
the exponentially scaling many-electron problem that is 
suitable for describing strong electron correlation. The 
2-RDM methodology has been applied in quantum chem- 
istry and condensed-matter physics to studyin g n rany- 
electron molecules and their reactions [2l|, l23ll33. 35 1, 
quantum phase transitions f36l , [37| , quantum dots 38 1 , 
molecular clusters )39 ^ 40| and spin systems like the Hub- 
bard [4l| and Ising [42| models. While new wavefunction 
methods for strong correlation are being actively devel- 
oped ;43.-.50i] ■ traditional wavefunction methods are lim- 
ited to linear combinations of approximately 10^ molec- 
ular configurations. 



2 



In this article we present a new approach to improv- 
ing the accuracy of energies from the 2-RDM method 
through rank restriction of the iV-representabihty con- 
ditions. The accuracy of the 2-RDM calculations with 

2- positivity can be dramatically enhanced through the 
addition of 3-positivity conditions [13, [13, Ell, H^, but 
these conditions significantly increase the total compu- 
tational cost of the calculations. Rather than turning to 

3- positivity, we propose to improve 2-positivity more di- 
rectly without an increase in its computational cost. The 
central idea is that approximate N -representability condi- 
tions like the 2-positivity constraints should be combined 
with less flexibility in the 2-RDM than more stringent 
N -representability conditions like 3-positivity. 

One systematic approach to controlling the flexibility 
of the 2-RDM is to restrict its rank or the rank of one of 
its metric matrices. We can motivate the selection of the 
rank by examining the ranks associated with model wave- 
functions including mean- field i53j and antisymmetrized 
geminal power (AGP) [3, [H-ilwavefunctions. Impor- 
tantly, because the 2-positivity conditions yield the exact 
ground-state energy of any AGP Hamiltonian — that is, a 
Hamiltonian whose ground state is described by an AGP 
wavefunction 

the AGP 2-RDM provides a lower 
bound on the optimal rank of the particle-hole form of the 
2-RDM. The resulting rank-restricted iV-representability 
conditions yield significantly improved ground-state en- 
ergies at a slightly lower computational cost than unre- 
stricted 2-positivity conditions. 

After theoretical results are presented in section [III il- 
lustrative applications are made in section [TlT] to comput- 
ing ground-state energies for a set of molecules in several 
basis sets as well as bond stretching of hydrogen fluo- 
ride and diatomic nitrogen. Section IIVI provides a brief 
discussion and concluding remarks. 

II. THEORY 

The energy is expressed as a functional of the 2-RDM 
in sections III Al and the iV-representability constraints, 
known as 2-positivity conditions Q , are reviewed in 
sections III Bl In section III C II we derive the maximum 
rank of the particle- hole matrix for two model wave- 
functions, the Hartree-Fock wavefunction and the AGP 
wavefunction. For the AGP wavefunction the maximum 
rank for each block of the spin-adapted particle-hole 
matrix is also derived. Finally, in section lTl C 2l we extend 
a large-scale algorithm for SDP [l^, [2^ to support rank 
restriction. 



A. Energy functional 

Because electrons are indistinguishable with pairwise 
interactions, the energy of any iV-electron quantum sys- 
tem can be expressed as a linear functional of the two- 
electron reduced Hamiltonian matrix and the two- 



electron reduced density matrix (2-RDM) [ll-Q 

E=j2 'Ki:rDi^:i (1) 

E = Ti{^K ^D), (2) 

where the indices denote spin orbitals 0p in a finite one- 
electron basis set {4>p}. The elements of the two-electron 
reduced Hamiltonian matrix are 

^K^J^-L-^KPS! + 'V!'f, (3) 

where matrices and contain the one- and two- 
electron integrals respectively, and the elements of the 
2-RDM are 

'D'^J = {^\alala,a,\^) (4) 

where at (a^) is a creation (annihilation) operator in sec- 
ond quantization that creates (annihilates) an electron 
in spin orbital cjjp and 5" represents the A'^-electron wave- 
function. 



B. A'^-representability conditions 

Because not every two-electron density matrix is rep- 
resentable by an iV-electron density matrix, the 2- 
RDM must be constrained by iV-representability condi- 
tions 2, 4-6l.l8l-l23|. A systematic hierarchy of constraints 
is furnished by the p-positivity conditions [1, H, [111, [O] ■ 
The 1-positivity conditions, constraining the one-particle 
and the one-hole RDMs to be positive semidef- 
inite, corresponds to restricting the eigenvalues of the 
1-RDM rip, known as natural occupation numbers, to lie 
between zero and one n„ g [0, 1], which enforces the Pauli 
principle. Coleman first proved that these relatively 
simple conditions plus the usual trace, Hermiticity, and 
antisymmetry constraints in the definition of a density 
matrix are not only necessary but also sufficient for the 
1-RDM to represent an TV-electron density matrix. 

The 2-positivity conditions 0, Q , providing necessary 
constraints on the 2-RDM, constrain the following three 
metric matrices to be positive semidefinite; 

'^DtO (5) 
^QhO (6) 
^G h 0, (7) 

where the metric matrices ^D, ^Q, and ^G correspond to 
the probability distributions for two particles, two holes, 
and one particle and one hole. In second quantization 
the elements of these matrices are expressible as 

'x^:? = {^\^Cp,,''cl\^), (8) 

where 

^Cp,, = alal (9) 

QCp.g = OpCLq (10) 

^Cp.q = alcLq. (11) 
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All three metric matrices contain equivalent information 
in the sense that rearranging the creation and annihi- 
lation operators produces linear mappings between the 
elements of the three matrices [l|, [2, [ll| ; particularly, the 
2-hole RDM and the particle- hole RDM can be 
written in terms of the 2-particlc RDM ^I? as follows 



and 



2 TP- q 



'2QP,q _ 1 



2jjP,q 



2DP,t 



(12) 



(13) 



where and are the one- and two-particle iden- 
tity matrices and A denotes the Grassmann wedge prod- 
uct [1^ [13] ■ While all three matrices are interconvert- 
ible, the nonnegativity of the eigenvalues of one matrix 
does not imply the nonnegativity of the eigenvalues of 
the other matrices, and hence, each semidefinite con- 
straint in Eqs. ([S]), ([HI), and ([7]) provides an important 
A^-representability condition. 



C. Rank restriction 

1. Model wavef unctions 

The best known model wavefunction is the mean- 
field wavefunction introduced by Hartree, Fock, and 
Slater [s^ . In first quantization the A^-electron Hartree- 
Fock wavefunction can be expressed as 

«'HF = 0l(l)A02(2)A...A(/.Ar(iV), (14) 

while in second quantization it can be written as 



(15) 



where |0) denotes the vacuum state, the state without 
any electrons. The rank of the particle-hole 2-RDM (or 
^G), whose elements are given in Eq. equals the num- 
ber of linearly independent A'^-electron functions \fij) 
having the form 



(16) 



For the Hartree- Fock wavefunction the set {\fi.j)} con- 
tains the wavefunction IvI/hf) itself as well as (r — N)N 
functions from all single excitations of I^'hf)- Hence, 
the rank of ^G from a Hartree- Fock wavefunction is 
(r — N)N + 1. For a wavefunction to describe a cor- 
related A^-clectron system in r spin orbitals its particle- 
hole 2-RDM ^G must have a rank strictly larger than 
(r- iV)iV + 1. 

A flexible model wavefunction with electron correla- 
tion is the antisymmetrized geminal power (AGP) wave- 
function P. [53 - [63 |. also known as the projected Bardeen, 
Cooper, and Schrieffer (BCS) wavefunction, which can be 



employed to model Cooper pairing in superconductivity. 
The A^-electron AGP wavefunction in first quantization 
can be written as 

*AGP - <?(1, 2) A 5(3, 4) A ... A g{N - 1, iV), (17) 

where g(l,2) is a two- electron function (or geminal) in 
contrast to the set of one-electron orbitals {(f)i]- In sec- 
ond quantization we can define the AGP wavefunction as 
a projection of the BCS wavefunction onto the A^-electron 
space 



I^AGI 



P, 



TV 



r/2 



|0), (18) 



where |0) is the vacuum state and P^r is the projection 
operator that projects the BCS wavefunction onto the 
Hilbert space of iV-electron wavefunctions. A key fea- 
ture of the AGP wavefunction is the special pairing of 
orbitals [2j], which we denote by and —i for i S [1, r/2]. 
In the study of superconductivity this pairing is employed 
to model the observed Cooper pairing of the momenta of 
electrons. 

The rank of ^G, again equaling the number of lin- 
early independent iV-electron functions in Eq. (|16p . 
can be determined for AGP from the pairing of orbitals. 



For AGP the functions 
classes i57il: 



\fi,j) can be divided into two 



1/^) =PmI*agp) 

1/^) =Q.,,|*AGP), 



(19) 
(20) 



where the P^^j are projection operators and Qi j are oper- 
ators whose adjoint operators annihilate the AGP wave- 
function, that is 



Qt jvI/AGp) =0. 



(21) 



Specifically, when the 7^ are not more than doubly de- 
generate, r/2 linearly independent functions \f['^) arise 
from the projectors 



P ,; = a] a. 



(22) 



and r(r — 2)/2 linearly independent functions \ f^j) arise 
from the Qi^j operators whose adjoints are 



Q\.j = liO-l^] - sign(ij)7jaLja_ 



(23) 



where i,j G [— r/2, r/2]\[0] with i ^ j and i 7^ —j and 
sign(«j) returns the sign of the product of i and j. 

The fact that each of the r(r — 2)/2 operators Q\ j an- 
nihilates the AGP wavefunction follows from the pairing 
property of the orbitals 0, E^l- From the definition of 
the AGP wavefunction in Eq. ifTSj) , it can be seen that in 
each Slater determinant contributing to the AGP wave- 
function both orbitals in a pair, i.e. and (jj-ii ^'^^ 
either occupied or unoccupied. Furthermore, each pair 
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in the wavefunction is weighted by a corresponding ele- 
ment of the vector 7. Hence, the actions of the operators 
'jiolaj and ^ja}_-a^i on the AGP wavefunction are al- 
ways equal or opposite in sign depending on whether the 
function sign(i, j) is equal to or -1, which proves the 
result. The numbers of linearly independent and 

\f?j) ^'^^^ t)e less than their maximum values of r/2 and 
r(r — 2)/2 if the 7^ where 7^ = 7_i are more than doubly 
degenerate. Such a case occurs when the AGP wavefunc- 
tion reduces to the Hartree-Fock wavefunction and the 
numbers of linearly independent jj/'j) and \f^j) become 
1 and N{r ~ N), respectively. Gonsequently, for an AGP 
wavefunction the maximum rank of is r(r— 2)/2 + r/2 
or r{r — l)/2. 

In electronic calculations, when the expectation value 
of the z-component of the spin operator (Sz) vanishes, 
the basis functions of the metric matrix can be spin 
adapted to produce a block diagonal matrix with four 
blocks [65l|. The four blocks correspond to the following 
four ^Cp,q operators: 



V2 



(1:-1) 



— =(a4 



Gc^+^ = dl 

I. J la 



(24) 
(25) 
(26) 
(27) 



where the bar above the index refers to the spatial part 
of the orbital, the spin part of each orbital is denoted as 
either a (-(-1/2) or (3 (—1/2), and the upper right indices 
of C^-j denote the square of the total spin and the 
z-component of the total spin for the two-electron oper- 
ators. If the pairing within the AGP ansatz is taken to 
be between spin orbitals sharing the same spatial compo- 
nent, the AGP wavefunction in Eq. (|18|) can be re- written 
with ia and «/3 replacing and —i. 

To determine the rank of the ^G spin blocks, we can 
spin adapt the projection operators in Eq. (I22[) and the 
adjoint of the annihilation operators in Eq. (I23|) . respec- 
tively, to obtain 



and 



p(0,0) ^ J_^p_ _ ^P- - \ 



(28) 

(29) 
(30) 
(31) 
(32) 



All Vs spin-adapted projection operators contribute to 
the (0,0) spin block of ^G, and rs{rs — l)/2 Q-type op- 
erators contribute to each of the (0,0), (1,-1), (1,0), 



and (1,-1-1) spin blocks where the number of spatial 
orbitals equals one-half the number r of spin orbitals. 
Hence, the rank of the singlet spin block (0, 0) of ^G is 
rs{r s + l)/2, and the ranks of the three triplet spin blocks 
of ^G are rs(r^ - l)/2. When {S^} = 0, all three triplet 
blocks are identical ^5|. 



2. Semidefimte programming 

The variational 2-RDM method with 2-positivity con- 
ditions minimizes the ground-state energy as a 2-RDM 
functional 



minimize 



E{x) 



T 
C X 



(33) 



where the vector c contains information about the quan- 
tum system in the form of the two-electron reduced 
Hamiltonian in Eq. ^ Q and the vector x contains the 
three different metric-matrix forms of the 2-RDM whose 
elements are given in Eq. (|H]). Because the three metric 
matrices in x are interrelated by linear mappings 



Ax 

and constrained to be positive semidefinite 

M{x) = 




(34) 



(35) 



where the operator M maps the vector x to a matrix, 
the energy minimization constitutes a special type of 
constrained optimization known as semidefinite program- 
ming (SDP) iMl m-IH. SDP is a generalization of 
linear programming from linear scalar inequalities to lin- 
ear matrix inequalities. 

Second-order algorithms for SDP, developed in the 
1990s [30, 31], have an expensive r^^ scaling [l^, [HI in 
floating-point operations when applied to variational 2- 
RDM calculations with 2-positivity constraints. Zhao et 
al. [12| introduced a dual formulation of the 2-RDM op- 



timization that decreased the computational scaling to 
r^^, and one of the authors developed two first-order al- 
gorithms, a matrix- factorization method and 
a boundary-point method '2^, that reduce the floating- 
point operations to and the memory requirements from 
to r^. Gances, Stoltz, and Lewin [15], who studied a 
dual formulation of the SDP problem, confirmed the effi- 
ciency of the matrix factorization method, and Verstichel 
et al. [22| introduced a first-order algorithm, based on 
interior-point methods. 

For the rank-restricted A^-representability conditions 
the SDP optimization must be modified to include rank 
restriction of the particle- hole ^G metric matrix within 
M . In the matrix-factorization method the solution ma- 
trix M is explicitly constrained to be p ositive semidefinite 
by a matrix factorization [ij, [3 126| : 



M = RR*. 



(36) 
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Importantly, the rank of M or any of its subblocks can 
be readily constrained to an integer q by restricting the 
number of columns of i? to g where q is less than the 
dimension of the square matrix M . With this flexibility 
we can solve SDP problems both with and without rank 
restriction. If the rank of a block in R is restricted to 
an unphysical value such as an integer less than the rank 
correspond to a Hartree-Fock model wavefunction, the 
algorithm generally will not converge. Otherwise, con- 
vergence of the rank-restricted SDP is similar to that of 
the unrestricted SDP. 



III. APPLICATIONS 



where Ns = N/2. 

(4) Contraction of the spin-adapted 2-RDM [65] onto 
the 1-RDM: 



1) 
1) 



(14) 

i.k\j.k 
(0,0) 



(40) 



(5) The 2-positivity conditions [Eqs. (l5]-[7])], on three 
different representations of the 2-RDM whose elements, 
given in Eq. are related by the linear mappings in 
Eqs. ^ and (HH). 

(6) In the case of rank restriction, the rank of the 
particle-hole matrix ^G^'^'"^ is restricted. 



After an overview of computational details and a sum- 
mary of A^-representability conditions, we present re- 
sults of the rank-restricted variational 2-RDM method 
for molecules at both equilibrium and non-equilibrium 
geometries. 

A. Computational details 

The variational 2-RDM method with 2-positivity and 
rank-restricted 2-positivity conditions is illustrated with 
calculations on several molecules at equilibrium and non- 
equilibrium geometries in minimal Slater-type orbital 
(ST0-6G) W], double-zeta (DZ) [g^], and correlation- 
consistent polarized double-zeta (cc-pVDZ) [6^ basis 
sets. Non-equilibrium geometries are obtained from the 
Handbook of Chemistry and Physics ^69^] , all core orbitals 
are double occupied (frozen), and the molecules are in 
singlet states. The calculation of one- and two-electron 
integrals and full configuration interaction (FCI) is im- 
plemented in the quantum chemistry package GAMESS 
(USA) 0. 

B. Summary of AT-representability conditions 

Variational RDM ground-state energies are computed 
with the first-order nonlinear SDP algorithm devel- 
oped by Mazziotti [11 [H, i^. The following N- 
representability conditions are enforced: 

(1) Hermiticity of the 2-RDM: 

= 'D^i • (37) 

(2) Antisymmetry of upper and lower indices 

= = - 'Dli = 'Di:l , (38) 

is enforced by antisymmetrized basis functions = 
l/\/2((?!)ij - (pj^i). 

(3) Trace conditions on the spin-adapted blocks of the 
2-RDM 

T,( 2^(1,0)) =N,{N,-1) 

Tr( 2^(0,0)) ^Ns{Ns + l) , ^ ' 



C. Results 

Two sets of ]V-representability constraints are imposed 
in the calculations shown in Tables Hl llIII and Figs. [T] and [2] 
(j) 2-positivity conditions without rank restriction, la- 
beled full rank, (ii) 2-positivity conditions plus rank re- 
striction, labeled theoretical rank, in which the rank of 
the 2(^(0,0) yock of the particle-hole matrix is restricted 
to rs{rs + l)/2 — its maximum value from a model AGP 
wavefunction. 

For a variety of molecules and basis sets Table H] shows 
the percentage of the correlation energy recovered by 
the variational 2-RDM method with full-rank and rank- 
restricted iV-representability conditions. Rank restric- 
tion significantly improves the percentage of correlation 
energy for all molecules and basis sets. For CO the 
2-positivity conditions without restriction yield 108.5%, 
119.3%, and 121.7% of the correlation energy in ST0-6G, 
DZ, and cc-pVDZ basis sets while these condition with 
rank restriction yield 101.2%, 102.5%, and 101.1% of the 
correlation energy. Even though the theoretical rank in- 
creases dramatically with basis-set size, the percentage 
of the correlation energy recovered remains nearly con- 
stant. Because the rank restriction limits the flexibility 
of the 2-RDM, we observe that the computed energies 
with rank restriction are neither consistently above nor 
below the FCI energy. 

Dissociation of the triple-bonded nitrogen molecule N2 
provides a classic case of strong electron correlation. Ta- 
ble [n] and Figure [T] present the potential energy curve 
of N2 in the cc-pVDZ basis set from the variational 2- 
RDM method with and without rank restriction. At 
R = 1.485 A, while the 2-RDM method with the fuh rank 
recovers 119.6% of the correlation energy, the 2-RDM 
method with the theoretical rank yields 99.6% of the cor- 
relation energy. At this distance in a region of the poten- 
tial energy curve where the spins are recoupling, some- 
times known as the spin recoupling region, the error from 
the rank-restricted 2-RDM method is only +0.002 a.u. 
relative to FCI. Figure [T] shows that the potential energy 
curve from the rank-restricted 2-RDM method closely 
agrees with the curve from FCI in a large region sur- 
rounding the equilibrium geometry. The largest errors 



6 



TABLE I. The percentage of the correlation energy recovered by the variational 2-RDM method with full-rank and rank- 
restricted Af-representability conditions is shown for the molecules CO, N2, H2O, HP, and NO^ in a variety of basis sets. 







Pull 


Theoretical 


Pull CI 


Correlation 


% Correlation Energy 


Molecule 


Basis Set 


Rank 


Rank 


Energy (a.u.) 


Energy (a.u.) 


Pull Rank 


Theoretical Rank 


CO 


STO-6G 


64 


36 


-112.443174 


-0.139676 


108.5 


101.2 




DZ 


256 


136 


-112.893590 


-0.208672 


119.3 


102.5 




cc-pVDZ 


676 


351 


-113.054884 


-0.305767 


121.7 


101.1 


N2 


STO-6G 


64 


36 


-108.699813 


-0.158189 


107.6 


97.7 




DZ 


256 


136 


-109.104172 


-0.226029 


119.1 


102.8 




cc-pVDZ 


676 


351 


-109.278339 


-0.329000 


121.5 


100.0 


H2O 


STO-6G 


36 


21 


-75.728838 


-0.050041 


104.0 


101.7 




DZ 


144 


78 


-76.141153 


-0.132021 


111.8 


102.9 




cc-pVDZ 


529 


276 


-76.241677 


-0.214915 


116.7 


107.3 


HP 


ST0-6G 


36 


21 


-99.526353 


-0.026196 


100.0 


99.7 




DZ 


100 


55 


-100.146049 


-0.124147 


107.2 


100.5 




cc-pVDZ 


324 


171 


-100.228652 


-0.209363 


115.2 


97.8 


N0+ 


ST0-6G 


64 


36 


-128.637594 


-0.241971 


108.7 


102.4 




DZ 


256 


136 


-129.060275 


-0.315068 


115.3 


103.6 



TABLE II. The percentage of correlation energy along the potential energy curve of the nitrogen molecule N2 in the cc-pVDZ 
basis set is reported from the variational 2-RDM method with and without rank restriction. At R = 1.485 A in a region of 
the potential energy curve where the spins are recoupling, sometimes known as the spin recoupling region, the error from the 
rank-restricted 2-RDM method is only -1-0.002 a.u. relative to PCI. 



Bond 


Pull CI 


Correlation 


% Correlation 


Energy 


Length (A) 


Energy (a.u.) 


Energy (a.u.) 


Pull Rank 


Theoretical Rank 


0.80 


-108.664476 


-0.257118 


120.7 


97.5 


1.1208 


-109.282139 


-0.332762 


120.1 


98.7 


1.175 


-109.275424 


-0.347604 


120.1 


98.7 


1.475 


-109.141160 


-0.442743 


119.6 


99.6 


1.85 


-109.008801 


-0.590167 


116.4 


97.5 


2.225 


-108.970662 


-0.748815 


111.3 


96.2 


2.6 


-108.963937 


-0.874616 


108.2 


95.5 


2.975 


-108.962249 


-0.963251 


106.8 


95.4 



from rank restriction occur at significantly stretched ge- 
ometries where strong spin entanglement increases the 
actual rank of the particle- hole matrix. In contrast, 
as observed in previous work, the 2-RDM method with- 
out rank restriction has its largest errors in the spin- 
recoupling region of the potential energy curve. One mea- 
sure for the potential curve's shape is the non-parallelity 
error, the difference between the largest error and the 
smaller error along the curve relative to FCI. While the 
2-RDM methods with and without rank restriction have 
similar non-parallelity errors over the whole curve shown 
in Fig. [U in the region R G [0.8, 1.85] the rank restric- 
tion improves the non-parallelity error from 0.044 a.u. to 
0.013 a.u. 

Due to the high electronegativity of fluorine, the disso- 
ciation of the hydrogen fluoride molecule illustrates the 
breaking of a polar covalent single bond. Table IIIII and 
Figure [2] present the potential energy curve of HF in the 
cc-pVDZ basis set from the variational 2-RDM method 
with and without rank restriction. At all bond lengths 
the rank restriction significantly reduces the error in the 
percentage of the correlation energy relative to FCI. At 



R = 1.95 in the spin- recoupling region the rank restric- 
tion reduces the error from —44.5 a.u. to —3.1 a.u. In 
the region R G [0.7, 2.8] shown in Fig. [2] the rank restric- 
tion decreases the non-parallelity error from 0.027 a.u. to 
0.024 a.u., and in the region R G [0.7,1.95] it decreases 
the non-parallelity error from 0.014 a.u. to 0.010 a.u. 
In contrast to the triple-bonded N2, the singly bonded 
hydrogen fluoride at highly stretched geometries has en- 
ergies from the rank restriction below those from FCI. 



IV. DISCUSSION AND CONCLUSIONS 

Variational minimization of the ground-state en- 
ergy as a function of the 2-RDM [9l-[23]. con- 
strained by A^-representability conditions, provides a 
polynomial-scaling approach to studying strongly cor- 
related molecules without computing the many-electron 
wavefunction. Here we have introduced a new approach 
to enhancing necessary conditions for 7V-representability 
through rank restriction of the 2-RDM. Applications 
were made to molecules at both equilibrium and non- 
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TABLE III. The percentage of correlation energy along the potential energy curve of hydrogen fluoride in the cc-pVDZ basis 
set is reported from the variational 2-RDM method with and without rank restriction. At R — 1.95 A in the spin-recoupling 
region the rank restriction reduces the error from -44.5 a.u. to -3.1 a.u. 



Bond 


Full CI 


Correlation 


% Correlation 


EnerETv 


Length (A) 


Energy (a.u.) 


Energy (a.u.) 


Full Rank 


Theoretical Rank 


0.70 


-100.129860 


-0.199411 


115.3 


96.1 


0.9161 


-100.228633 


-0.209189 


115.2 


97.6 


1.2 


-100.181953 


-0.222427 


115.3 


98.0 


1.3 


-100.157836 


-0.227830 


115.3 


98.3 


1.5 


-100.113798 


-0.240945 


115.5 


97.2 


1.95 


-100.052917 


-0.283544 


115.7 


101.1 


2.8 


-100.026420 


-0.369590 


115.6 


104.7 



-108 








-108.2 








-108.4 






+ RHF 

Theoretical Rank 


-108.6 






- FCI 

□ Full Rank 


-108.8 
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1 1.5 2 2.5 



N-N Distance 

FIG. 1. The potential energy curve of the nitrogen molecule 
N2 in the cc-pVDZ basis set from the variational 2-RDM 
method with and without rank restriction. Results are com- 
pare with those from Hartree-Fock and FCI. 
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FIG. 2. The potential energy curve of the hydrogen fluoride 
molecule in the cc-pVDZ basis set from the variational 2- 
RDM method with and without rank restriction. Results are 
compare with those from Hartree-Fock and FCI. 



equilibrium geometries. 

An important set of 7V-representability conditions on 
the 2-RDM is the 2-positivity conditions, which restrict 
the probability distributions of two particles (^D), two 
holes (^Q), and a particle-hole pair i^G) to be nonnega- 
tive. In section III C II we derived the maximum rank of 
the particle- hole metric matrix for two model wave- 
functions, the Hartree-Fock and the AGP wavefunctions. 
The Hartree-Fock wave functions are a small subset of the 
AGP wavefunctions, and hence, their particle-hole matri- 
ces have a maximum rank (r — N)N + 1 which is strictly 
less than the maximum rank r{r — l)/2 of the AGP 
particle-hole matrices. Because the 2-positivity condi- 
tions constrain AGP Hamiltonians — that is, Hamiltoni- 
ans with AGP ground-state wavefunctions — to yield the 
exact ground-state energies, the rank of the AGP ma- 
trix provides a minimum rank for the molecular particle- 



hole ^G matrix within variational 2-RDM calculations of 
general systems. Selecting a smaller rank for ^G would 
render the variational 2-RDM method inexact for AGP 
Hamiltonians. Unlike the case in Hartree-Fock theory, 
the maximum rank of within the AGP ansatz is inde- 
pendent of the number N of particles, which reflects its 
independence from a reference determinant wavefunction 
and hence, its ability to treat strong electron correlation. 

The variational 2-RDM method with rank-restricted 2- 
positivity conditions was applied to computing the ener- 
gies and 2-RDMs for a variety of molecules at equilibrium 
geometries as well the potential energy curves of the ni- 
trogen and hydrogen fluoride molecules. Specifically, the 
rank of the singlet spin block of the particle-hole ma- 
trix was restricted to its maximum value from an AGP 
wavefunction rs{rs + l)/2 with = r/2. The rank- 
restricted conditions were implemented through an ex- 
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tension of the first-order matrix-factorization algorithm 
for large-scale SDP. The results demonstrate that rank 
restriction significantly improves the accuracy of com- 
puted energies. For example, the percentages of corre- 
lation energies recovered for HF, CO, and N2 improve 
from 115.2%, 121.7%, and 121.5% without rank restric- 
tion to 97.8%, 101.1%, and 100.0% with rank restric- 
tion, respectively. The improvement occurs at equilib- 
rium and non-equilibrium geometries and across basis 
sets. Computationally, the rank-restricted conditions are 
slightly less expensive than the full 2-positivity condi- 
tions. Rank restriction removes degrees of freedom that 
are not sufhciently constrained by the 2-positivity con- 
ditions without sacrificing the method's ability to treat 
strong electron correlation, as seen in the bond dissocia- 



tion of N2 . Although further research is needed to study 
the method in larger molecules such as polyaromatic hy- 
drocarbons [U 111] and firefly luciferin [l^l , the present 
results indicate that rank restriction is a promising ap- 
proach to improving the 2-positivity conditions within 
the variational 2-RDM method without increasing com- 
putational cost. 
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